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Abstract 

By using photon pairs created in parametric down conversion, we report on an experiment, which 
demonstrates that measurement can recover the quantum entanglement of two qubit system in a 
pure dephasing environment. The concurrence of the final state with and without measurement 
are compared and analyzed. Furthermore, we verify that recovered states can still violate Bell's 
inequality, that is, to say, such recovered states exhibit nonlocality. In the context of quantum 
entanglement, sudden death and rebirth provide clear evidence, which verifies that entanglement 
dynamics of the system is sensitive not only to its environment, but also on its initial state. 
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Quantum entanglement, as a unique feature without a classical counterpart of many-body 
system, has been instrumental in studying fundamental aspects of quantum physics Q as well 
as being central in practical applications in the areas of quantum computation and quantum 
cryptography j^lj]. Today, sources of entangled states can be prepared in various kinds 
of physical systems Entanglement, even within multi-particle and multi-dimensional 
systems [7|,|8j, can be implemented, although, these are flimsy and are subject to unavoidable 
degradation, which is caused by interactions with their environments [9_|, |l0|. To have any 
practical value in quantum computation and communication, long distance nonlocality and 



extended coherence storage and rebirth 



ll| of entangled states have become important 



focal points of research around the world. The key issue behind solving these problems is 
in determining the dynamical behavior of entanglement within the system's environment, 
something, which to date, has not been well understood. Commendably, a factorization 
law, which describes the entanglement dynamics under a one-side noisy channel has recently 
)een proposed [3] and has subsequently been verified in two independent experiments 
131 1 . Moreover, it has been discovered that entanglement evolution is not only related to 
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environmental factors but is also sensitive to its initial state |14J . Il5j . 

Quantum measurement, that feature, which distinguishes the quantum from the classi- 
ca. regimes H, is often interpreted within the orthogonal projection mode, given by Von 
Neumann [18j, but has been reexpressed in the past 30 years or so, more and more within 
the framework of quantum decoherence theory 19] (for reviews see jsl, In that setting, 
a complete quantum measurement is divided into two stages: The first stage corresponds to 
the entanglement of the information-carrying qubit of the measured quantum system with 
the record bit of the measurement apparatus; the second stage corresponds to the deco- 
rerence of the detector-system combination, which occurs in an uncontrollable environment 
17j. The latter will covert the density matrix that describes the combination into diagonal 
form, which physically represents the spreadings of the quantum information contained in 
the combined system into the uncontrollable environment and essentially turns the quan- 
tum measurement into a classical one. Because of the uncontro liability and unavoidability of 
environment-induced decoherence that features prominently in this second stage, it is impos- 
sible to recover quantum information once the quantum measurement has been completed. 
For example, if a conventional experiment to measure the polarization of a single photon has 
completed, to retrieve that qubit of quantum information is as difficult as extracting it from 
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the single photon detector and the observer's brain. Fortunately, quantum measurement 
need not be so catastrophic and can be implemented step by step, even partially [20]. It 
has been pointed out that, within the first stage of quantum measurement, recovery can be 
brought about by the coherence of a single qubit that has dissipated into a non-Markovian 
environment 2l| . 



On the basis of this work by Xu et al Xu we experimentally prove that subsequent 



quantum measurement, which erases 



16] the path information introduced in the previous 



measurement can recover the entanglement that has been degraded in a non-Markovian envi- 
ronment and that even rebirth of entanglement after entanglement sudden death (ESD) 11] 
can be effected. In this paper, we describe how we can change the entanglement evolution by 
using some specified operations on the state before or during the interaction, more precisely, 
two sequential measurements can recover the entanglement and, further more, preserve it. 

The experimental setup is shown schematically in Figure 1. Two 0.5 mm thick beta- 
barium-borate (BBO) crystals, cut at 29.18° for type-I phase matching and aligned so their 
optical axes are perpendicular to each other, are pumped by using focused ultraviolet (UV) 
pulses polarized at 45°, which are frequency doubled from a Ti:sapphire laser with the cen- 
ter wavelength mode locked at 800 nm (with 130 fs pulse width and a 76 MHz repetition 
rate). Degenerate polarization-entangled photon pairs at 800 nm are generated by sponta- 
neous parametric down conversion (SPDC) at a 3° angle with the pump beam 22|. By 
compensating the time difference between H- and ^/-polarized components with birefrin- 
gent elements (LiNbOs and YVO4), one of the maximal polarization-entangled states, the 
well-known Bell state 30], can be produced with high fidelity. This initial state can be 
mathematically written as 

\<f>} = ^(\HH} + \VV}), (1) 

where H and V represent the horizontal and vertical polarizations, respectively, while the 
two elements in the Dirac ket label the photons states, the left in path a and the right in 
path b. 



lab 



e birefringent elements, 



241 ] . In our experiment, 



Decoherence due to the environment is simulated by contra, 
which can couple the photon's frequency with its polarization 23, 
we use quartz plates Ql with thickness L\ and Q2 with thickness L2, to simulate this 
decoherence aspect. The optical axes of the plates are both horizontally set. 
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FIG. 1: (color online) Scheme of the experimental setup. DE, decoherence evolution denoted by 
a dashed pane; Measurement apparatus (M) is denoted by two gray boxes. PBS, polarizing beam 
splitter; LI and L2, lens; PA, polarization analyzer; SPAD, single photon detector; BD1 and BD2, 
beam displacing prism; The solid pane E is inserted in path a to prepare the partially entangled 
state. The parameter of these elements are provided in the text. 

At the beginning, we consider the evolution of the state given by equation (1) in such an 
environment, which assume, for simplicity, that the photon's frequency distribution is a 5 
function. The final state of the two photons after the photon in mode b passes through Ql 
and Q2 takes the following form 

\<j )l ) = -L(\HH) + e^\VV)), (2) 

where the parameter a is proportional to L{= L\ + L2) and Ub represents the frequency 
of the photon in path b. In our experiment, a = LAn/c where c is the velocity of light 
in vacuo and An = n — n e represents the difference between refractive indices of ordinary 
(n Q ) and extraordinary (n e ) light. Therefore, a deterministic relative phase is introduced 
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between \HH) and \VV) for a sing 



cording to the decoherence mode 



23 



e frequency, which differs for different frequencies. Ac- 



24J , the environment is actually composed of a photon's 
frequencies that are coupled to the information carriers (viz. polarization of photons) by 
means of birefringent elements. For photons with a frequency distribution f(uj), the overall 
state should equal the integral over the frequency distribution, which forms a less correlated 
state, which essentially destroys the coherence of the qubits [23j. There: 
in Eq. (2) should be replaced by the following reduced density operator 24] 



ore, the final state 



Px = 1 -{\HH)(HH\ + \VV){VV\ + k* b \HH)(VV\ + k b \VV)(HH\), (3) 

where the photon's frequency distribution in path b is normalized as / f(ub)dub = 1 and 
the nondiagonal coefficient fc& = J f(ub)e tauJb dub represents the decoherence parameter. 

Actually, for a decoherence environment composed of photon's frequency, the decoherence 
parameter kb is related to the frequency distribution function. In our experiment, this is 
taken to be a Gaussian function, that is, f(uj b ) = -j=^ exp ( — 4( - aJb ~^ } ° ) ) ; which is determined 
by the interference filters (IF) placed in front of the single photon detectors. The parameter 
Co>o is the central frequency, and a is the bandwidth. By working out the integral, we obtain 
k b = exp(— a 2 cr 2 /16 + iauo). 

For a quantitative analysis of the entanglement evolution in the experiment, a parameter, 



which represents the degree of entanglement, should be introduced. The concurrence 
which is widely used in studying two qubit states, is defined as 



26|, 



C(p) = max{0, y/Xx - ^ X 2 - V A 3 - VM, (4) 

where p is the density matrix of a two-qubit state in the canonical basis 
{\HH), \HV), \VH), \VV)}; Xi{i = 0, ...,4) are the eigenvalues in decreasing order of 
the Hermitian matrix p(a y <g> a y )p*{a y <g> a y ) with p*, which corresponds to the com- 
plex conjugate of p According to Eq. (3), for an initial Bell state input, we obtain 
C(px) = \kb\ = exp(— a 2 cr 2 /8). Thus, concurrence degrades exponentially and approaches 
zero as Ql and Q2 become thicker, which means the final state, after sufficient interaction 
time, evolves into the maximally mixed state without any remaining entanglement [solid 
line in Fig. 2(A)]. However, for some other frequency distribution, some unusual phenomena 
will arise, for example, entanglement collapse and revival 25]. 

Then, we consider the case with a measurement apparatus (M), which comprises two 
beam displacing prisms (BD1 and BD2) with horizontally-set optical axes, half- wave plate 1 



(HWP1) with optical axes set at 22.5°, which implements the Hadamard operation, HWP2 
with optical axes set at -22.5°, which implements \H) —¥ |+),|V) — > — \— ) (where, |+) = 
(| if) + \V))/V2, |-) = (|ff) - \V))/y/2), and HWP3 with perpendicularly set optical axes 
to implement the bit-flip operation. BD1 measures photon's polarization in H/V basis and 
introduces the path information as a probe bit. Before this measurement is completed, the 
second decoherence environment is inserted, and then, the path information is erased by 
BD2, which realizes a postselection of the recovered state. The output state after passing 
through Ql, BD1, HWP1 and Q2, can be written as 

(^\H)(\H) + e ia ^\V)))j + (±e ia ^\V)(\H) - e™|y))) n (5) 

where a,\ = LiAn/c and 0:2 — L 2 An/c, subscripts I and II denote the upper and lower paths 
between the two BDs, respectively. By subsequently erasing path information introduced 
in the previous measuring operation by HWP2, BD2 and HWP3, we obtain the final state 
represented as 

|02) = ^(1 + e ia ^ B )(\HH) + e ta ^ B \VV)). (6) 
Similar to the above treatment, the reduced density matrix of this final state is written as 
p 2 = \{\HH)(HH\ + \VV)(VV\ + k' b *\HH)(VV\ + k' b \W)(HH\), (7) 
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where 



,/ exp(iaiwo) r r , \2 2 , ■ 1 

K = -7— -, r -, 2-^— — rjexp -(«! + a 2 ) (7 + ia 2 u \ 

2[1 + cos(a 2 ^o) exp(— a2<77 16)J 

+ exp[— («! — a 2 ) 2 cr 2 — 10:2^0] + 2 exp[— a\a 2 

According to Eq. (4), the concurrence is C(p 2 ) = \k' b \. From the complicated form of k' b , here, 
the entanglement evolution is not as simple as that in the previous case, and it sensitive 
to the phase. Numerical analysis shows there will be entanglement recovery with increases 
in L2, as L\ remains fixed. The concurrence oscillates if L 2 is sufficiently thin, while the 
amplitude narrows to zero if L2 is thick enough. Because the maximal recovery point lies 
within the envelope composed of the integral L2 and the zero point of the phase is extremely 
hard to determine, here, we only consider the integral length of the quartz plates [solid line in 
Fig. 2(B) and Fig. 2(C)]. Although this is sufficient for studying entanglement recovery, under 
this consideration, the probability of success decreases exponentially from 1 to 0.5. More 
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significantly, the entanglement of the final state will remain unchanged with L 2 increasing 
when L2 is large enough, and no matter how thick L\ is, that is, to say, no matter how 
less entanglement remains, the entanglement between two particles can be recovered to a 
maximal expectation value 0.5 at L 2 = L\. 

In fact, the dynamics of entanglement in bipartite quantum systems not only is sensitive to 
their environment, but also is sensitive to their initial state [l^LLs)]. In Ref 25|, entanglement 
collapse and revival occur when the input biphoton is prepared in the form of a Werner state 
j^sj with a spectrum discretized within a Gaussian envelope. There will be no entanglement 
revival if the spectrum takes the Gaussian form in that experiment. However, a revival of 
the same initial state with a Gaussian spectrum can occur by inserting M in this experiment. 
This is explained as follows. 

Applying a Hadamard operation on the photon state in mode a of the maximal entangled 
state |0) and by allowing it to pass though a dephasing environment at H/V bases, we get 
the state 



i 3 ) = ~(\HH) + \HV) + e iaaU)a \VH) - e iaaUJa \VV)). 



(8) 



By integrating ^3) (^3] over all frequencies of the photon in mode a, the partially entangled 
input state can be mathematically expressed in the following density matrix 
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(9) 



where k a = j g{ui a ) ex.\){ia a ui a )dui a is the decoherence parameter in path a. The photon in 
mode b then passes through the decoherence environment with a measuring apparatus; the 
final state in the single frequency case can be written as 



|03> 



4 



\\HH) + e iaiu)b \HV) + e iaaUla \VH} - e 



\VV)\. 



(10) 



The band width of the interference filters used in our experiment is so narrow that we 
can integrate \4>'z){(ftz\ over the frequencies of the photon in mode a and b separately 29]. 
Therefore, if the photon frequencies are considered to have Gaussian distribution profiles, 
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the density matrix of the final state is written as 



4 



1 



/ i k' b * K -k* a k'A 

\—k a k' b —k a —k' b 1 J 



(11) 



According to Eq. (4), we obtain a concurrence of C{pz) = max{0, (k a + k' b + k a k' b — 1)/2}. For 
nonnegative (k a + k b + k a k' b — l)/2 with different k' b , entanglement, sudden death and re-birth 
will occur if the decoherence time in Ql and Q2 are changed [solid line in Fig. 2(D) with 
just integral lengths of quartz plates]. No matter how thick Ql is, entanglement between 
the two particles is found to indicate a full rebirth to an identical maximal value at L 2 = L\. 
That is, to say, full recovery with analogous levels can be achieved regardless of the time 
duration taken for the ESD. 

The experimental results are shown in Fig. 2. The characteristics of the IFs used before 
the single photon detectors are of bandwidth 3 nm and coatings at 800 nm. We can treat 
An = 0.01 for small frequency distributions. The maximally entangled state is prepared 
with a concurrence of 0.962 ± 0.029. In Fig. 2(a), the concurrence degrades exponentially and 
gradually tends to zero, which obeys the half-life law. Because there is no phase sensitivity, 
experimental results (dots) agreed well with theory within the error range. In Fig. 2(b), M 
is inserted at the point Li = 195Ao where the concurrence is 0.262 ±0.024. The maximally 
recovered concurrence in the experiment is 0.609 ± 0.026 at L 2 = 195A , and the concurrence 
is unchanged within the error range when L 2 exceeds 683A . In Fig. 2(c), the point at which 
M is inserted is L 2 = 390Ao where the concurrence tends to zero, that is, to say, there 
are few entanglement at this point. Incidentally, the maximally recovered entanglement 
measured in the experiment is 0.518 ±0.025 at about L 2 = 780Ao, which agrees well with 
theoretical predictions. In Fig. 2(d), the partially entangled input state with concurrence 
0.704 ±0.019 is prepared by inserting an HWP with optical axes set at 22.5° and quartz 
plates of thickness 98A with horizontally-set optical axes in mode a. We insert M at 
L\ = 390 Ao at which ESD has occurred, and there is no entanglement; entanglement rebirth 
occurs at about L 2 = 585Ao and then the entanglement collapses at about L 2 = 975Ao 
again. The concurrence is corrected to zero according to Eq. (4) when its measured value is 
negative. Entanglement can be reborn at a maximal value 0.276 ±0.013 in the experiment. 
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FIG. 2: Major results (dots) for entanglement dynamics in the experiment. The solid lines are 
theoretical predictions of concurrence. Two error types are considered: the shot noise error in the 
measured coincidence counts and the uncertainty in the settings of the angles of the wave plates used 
to perform the tomography j^j ■ (a) represents the entanglement evolution of a maximal entangled 
state in a pure dephasing environment without M; (b) and (c) with M inserted at L\ = 195 and 
L\ = 390, respectively; (d) represents the entanglement evolution of the partially entangled state 
in a dephasing environment with M inserted at LI = 390. An = 800 nm. 

Nonlocality, as a particular characteristic of quantum mechanics, has changed the view- 
point and methods in understanding nature at its fundamental level. Since it can be studied 



by the well known Bell inequality 



301 ]. a quantum state with nonlocal correlations could be 



a very useful feature to exploit in future quantum technologies 3l|. In our experiment, the 
maximally recovered entangled state in Fig. 2(b) and Fig. 2(c) can be proven to be nonlocal 
by the more convenient Clauser-Horne-Shimony-Holt (CHSH) inequality [32|], S < 2 for any 
local realistic theory, where 



S = E{0 X , 6 2 ) + E(0 X , & 2 ) + E(6[, 9 2 ) - E{9[, 9' 2 ] 



(12) 



9 



with 

F(B n \ _ gfr, e 2 ) + c(e±, 9j) - c(e u 9j) - cjej, e 2 ) 

1 2) C(0 U 9 2 ) + C(6i, Oi) + C{9 U 6£) + C(6i, 9 2 ) ' 

Here, 9^(9^), i = 1,2 represent the linear polarization setting in path a and path b separately 
and 9f = #j + 90°, i — 1, 2. By calculating the maximal value of S from the measured density 
matrix of the maximally recovered state, we get (9± = —15°,^ = 21°, 9 2 = 86°,9' 2 = —52°) 
in Fig. 2(B) and (9 1 = -82°, 9[ = 66°, 9 2 = -4°,9' 2 = 28°) in Fig. 2(C). Accordingly, the 
measured values of S are 2.336 ±0.003 and 2.210 ±0.003, which violate the local realism 
limit 2 by over 104 and 64 standard deviations, respectively. 

The extraordinary phenomenon of entanglement recovery induced by a measurement, 
can be understood in the quantum framework of a partial measurement and reversal 20 ]. 
BD1 measures the photon's polarization in H/V basis. Partial measurement in H/V basis 
in the upper path and a reversal operation with the same strength {33 1 in the lower path 
is implemented by BD2. The photons in the two dark ports of BD2 are abandoned, this 
can be considered as a non-response of the detector to the coincidence detection. Because 
martial measurement and reversal with the same intensity can restore the initial entanglement 



20| and the phase difference introduced by L 2 only changes the strength, we believe that 



entanglement is preserved with the increasing L 2 while the recovery is an optical spin-echo 
effect introduced by HWP1. 

To summarize, we have report on an experiment, which shows that entanglement can 
be recovered by a process of measurement followed by quantum eraser in a non-Markovian 
environment. Simultaneously, the maximally recovered states can be verified to violate the 
CHSH inequality with high standard deviations, which confirms theirs quantum character. 
This result can be used to eliminate the influence of dephasing. Another encouraging as- 
pect is that, even if the state has been thoroughly disentangled, that is, to say, ESD has 
occurred, the entanglement can still be revived from non-Markovian environments regard- 
less of decoherence time durations. This can be used to implement controllable recovery of 
entanglement. Entanglement dynamics considering other coding and decoding protocol can 
be studied in future work. 
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